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Abstract — Analysis  of  intraseasonal  (10-100  days)  oscillations  in  the  equatorial  mesosphere  and  lower 
thermosphere  (MLT)  is  presented,  based  on  over  five  years  of  velocity  data  acquired  by  a  radar  system  at 
Christmas  Island  (2°N,  157  W),  in  the  central  Pacific.  Strong  peaks  in  the  zonal  winds  are  found  at  periods 
of  ~  60  days,  ~  35 — 40  days,  and  ~  22-25  days.  These  peaks,  as  well  as  the  mean  annual  variations  of  the 
activity  within  the  various  period  ranges,  are  similar  to  30-60  day  and  20-25  day  oscillations  that  occur  in 
the  equatorial  troposphere.  Weaker  (but  nonetheless  clear)  periodicities  are  also  found  in  the  meridional 
winds  at  ~  60  days  and  ~  35  days.  A  strong  quasi-60-day  variation  is  detected  in  gravity-wave  variances, 
with  much  weaker  signals  at  ~  40  days  and  ~  25  days.  Strong  variations  in  diurnal  tidal  amplitudes  are 
observed  with  periods  of  ~  60  days,  ~  40  days,  and  ~  25  days. 

These  observations  lead  us  to  propose  the  following  explanation  for  the  observed  intraseasonal  variability 
of  the  equatorial  MLT  region.  Intraseasonal  cycles  in  tropical  tropospheric  convection  produce  intra¬ 
seasonal  variations  in  the  intensity  of  gravity  waves  and  nonmigrating  diurnal  tides  impinging  upon  the 
mesosphere.  This  accounts  for  the  intraseasonal  peaks  we  observe  in  gravity-wave  and  tidal  activity.  This 
intraseaonally  modulated  wave  activity  induces  similar  periodicities  in  the  wave-induced  driving  of  the 
zonal  MLT  flow,  which  in  turn  forces  the  observed  intraseasonal  peaks  in  the  zonal  MLT  winds.  If  this 
explanation  is  valid,  these  observations  provide  an  unusually  clear  example  of  the  driving  of  MLT  flow 
patterns  by  waves  emanating  from  tropospheric  systems,  and  highlight  the  importance  of  convectively 
generated  waves  in  understanding  the  dynamics  of  the  equatorial  middle  atmosphere.  ©  1997  Elsevier 
Science  Ltd.  All  rights  reserved 


INTRODUCTION 

In  recent  years,  high-quality  long-term  observations 
of  the  equatorial  mesosphere  and  lower  thermosphere 
(MLT)  have  increased  our  knowledge  of  this  region 
considerably.  Equatorial  wind  measurements  from 
both  the  Upper  Atmosphere  Research  Satellite 
(UARS)  and  ground-based  radars  have  provided  us 
with  a  detailed  climatological  picture  of  the  mean 
circulation  (e.g.,  Lieberman  et  at.,  1993;  Palo  and 
Avery,  1993;  Garcia  et  al.,  1997),  equatorial  waves 
(e.g.,  Vincent,  1993;  Canziani  et  al.,  1995),  normal 
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modes  (Wu  et  al.,  1993;  Harris  and  Vincent,  1993; 
Palo  and  Avery,  1996),  tides  (e.g.,  McLandress  et  al., 
1996),  and  gravity  waves  (e.g.,  Isler  and  Fritts,  1996; 
Connor  and  Avery,  1996). 

On  analyzing  radar  winds  from  a  mid-Pacific  equa¬ 
torial  site,  Eckermann  and  Vincent  (1994)  presented 
preliminary  evidence  of  strong  intraseasonal  zonal- 
wind  oscillations  (periods  ~  10-100  days)  at  heights 
between  about  80-98  km.  They  noted  that  this  varia¬ 
bility  in  many  ways  resembled  intraseasonal  oscil¬ 
lations  encountered  in  the  equatorial  troposphere 
(Madden  and  Julian,  1994),  but  argued  that  the  large- 
scale  tropospheric  wave  disturbances  associated  with 
this  activity  could  not  propagate  to  mesospheric 
heights.  However,  since  this  intraseasonal  activity 
modulates  convective  activity  within  the  Indian  and 
western  Pacific  regions  (Hartmann  and  Gross,  1988; 
Dunkerton  and  Crum,  1995),  Eckermann  and  Vincent 
(1994)  hypothesized  that  a  similar  intraseasonal  vari¬ 
ation  in  the  intensity  of  convectively  generated  gravity 
waves  might  result.  If  these  waves  propagated  to 
mesospheric  heights  and  dissipated,  similar  intra- 
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Locations  of  Various  Observing  Sites 


Fig.  1.  Location  of  Christmas  Island,  as  well  as  some  other 
ground-based  stations  which  provided  data  which  we  discuss 
later. 


seasonal  variability  might  then  be  transferred  into  the 
equatorial  MLT  flow. 

The  study  of  Eckermann  and  Vincent  (1994), 
however,  was  preliminary.  Its  principal  aim  was  to 
present  a  persuasive  case  for  the  existence  of  these 
oscillations,  and  to  offer  some  possible  explanations 
of  their  cause.  Here,  we  investigate  the  intraseasonal 
activity  in  these  data  further,  extending  both  the 
breadth  of  the  analysis  and  the  length  of  data  analy¬ 
zed.  We  also  compute  gravity-wave  and  tidal  activities 
from  higher-resolution  wind  data  to  investigate 
whether  they  correlate  in  any  way  with  the  intra¬ 
seasonal  wind  variations,  which  might  help  to  assess 
the  gravity-wave-based  theory  for  their  generation 
suggested  by  Eckermann  and  Vincent  (1994). 

DATA 

The  wind  data  used  in  this  study  were  acquired  by 
a  spaced  antenna  MF  radar  stationed  at  Christmas 
Island  (2°N,  157°W),  in  the  equatorial  Pacific  (see  Fig. 
1).  The  system  measures  horizontal  wind  velocities 
every  2  min  at  2  km  height  intervals  in  the  range  60- 
98  km  (Vincent  and  Lesicar,  1991).  The  raw  data  were 
initially  screened  using  an  outlier-rejection  algorithm, 
then  diurnally  averaged  to  give  daily  wind  estimates  at 
each  height,  as  described  by  Eckermann  and  Vincent 
(1994).  Eckermann  and  Vincent  (1994)  analyzed 
slightly  less  than  3  yr  of  data,  obtained  between  23 
January,  1990  and  21  December,  1992  (1064days), 
whereas  here  we  incorporate  data  up  to  21  May,  1995 


(1944  days  in  all).  These  data  provide  us  with  nearly 
continuous  time  series  of  daily  zonal  and  meridional 
winds,  recorded  at  2  km  height  intervals  in  the  60- 
98  km  height  range,  which  we  use  to  study  intra¬ 
seasonal  activity.  Since  diurnal  data  coverage  below 
~  75-80  km  is  incomplete,  tidal  effects  may  not  be 
fully  removed  from  the  daily-averaged  data  at  these 
lower  heights,  and  so  we  focus  mostly  on  the  data 
above  80  km. 

To  study  gravity- wave  and  tidal  activity,  we  use 
higher-resolution  time  series  which  were  formed  by 
averaging  the  screened  raw  data  into  4  min  segments. 
These  data  extend  to  the  end  of  1994  only.  Gravity- 
wave  variances  were  computed  by  first  band-passing 
the  full  5  yr  time  series  using  a  boxcar  filter  spanning 
the  12min-4h  band,  with  missing  data  linearly  inter¬ 
polated.  Thereafter,  variances  were  computed  by  aver¬ 
aging  the  band-passed  data  over  2  day  blocks,  with 
interpolated  regions  excluded  from  the  analysis. 
Modeling  by  Rastogi  et  al.  (1996)  indicates  that  hori¬ 
zontal-velocity  variances  at  periods  <  1  h  can  be  arti¬ 
ficially  enhanced  by  systematic  instrumental  effects. 
However,  given  the  ‘red’  nature  of  wave  frequency 
spectra,  variances  computed  over  the  12  min^l  h  band 
are  influenced  more  by  the  longer  periods  ( >  1  h) 
which  are  less  affected  by  these  biases  (Rastogi  et  al., 
1996).  Tidal  structures  were  fitted  using  least-squares 
harmonic  fits  over  a  6  day  fitting  window  which  was 
successively  forwarded  every  2  days.  These  procedures 
yielded  time  series  of  both  gravity-wave  variance  and 
tidal  amplitudes  every  2  days  for  the  full  5  yr  of  data. 


SPECTRAL  ANALYSIS 

We  begin  by  spectrally  analyzing  the  velocity  time 
series  at  a  selected  height,  following  and  extending  a 
similar  preliminary  analysis  by  Eckermann  and  Vinc¬ 
ent  (1994).  Figure  2  profiles  time  series  of  zonal  and 
meridional  velocities  at  a  height  z  =  92  km.  The  top 
row  of  plots  shows  the  raw  daily-averaged  data,  while 
the  row  below  it  shows  the  same  data  after  smoothing 
with  a  20  day  running  average. 

Next,  power  spectral  densities  of  the  daily-averaged 
velocity  time  series  were  computed  using  numerical 
Fourier- transform  methods.  As  in  Eckermann  and 
Vincent  (1994),  small  gaps  in  the  time  series  due  to 
instrument  downtime  were  interpolated  prior  to  Four¬ 
ier  transformation  using  a  cubic  spline  under  tension. 
Resulting  power  spectral  densities  are  plotted  on  the 
third  row  of  Fig.  2  as  a  function  of  the  harmonics 
of  the  fundamental  frequency  ay  =  (1 944  days)  a 
corresponding  period  scale  in  days  is  shown  on  the 
upper  axis  for  reference. 


Winds  over  Christmas  Island  (2°N,  157°W) 
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Fig.  2.  Spectral  analysis  of  time  series  at  a  height  of  92  km  of  zonal  (left  column)  and  meridional  (right 
column)  velocities.  The  first  row  of  plots  shows  the  raw  daily-averaged  time  series,  which  are  then  smoothed 
with  a  20  day  running  average  and  plotted  on  the  next  row.  Power  spectra  of  the  daily  data  are  shown  on 
the  third  row.  A  Hanning  filter  function  is  also  shown  with  a  dotted  curve  on  these  plots.  It  is  used  to 
band-pass  fluctuations  at  ~  60  days.  The  resulting  band-passed  fluctuations  are  shown  on  the  bottom  row 

of  plots. 
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These  spectra  show  similar  features  to  those  com¬ 
puted  from  more-limited  data  by  Eckermann  and 
Vincent  (1994).  The  zonal- velocity  spectrum  is  domi¬ 
nated  by  a  large  off-scale  peak  at  ~  180  days  due  to 
the  mesosphere  and  lower-thermosphere  semiannual 
oscillation  (MLTSAO),  as  well  as  smaller  subsidiary 
peaks  suggestive  of  slight  annual  (e.g.,  Vincent,  1993) 
and  quasi-biennial  influences.  Additionally,  as  noted 
by  Eckermann  and  Vincent  (1994),  there  is  rich  struc¬ 
ture  in  the  spectrum  within  the  10-100  day  period 
range:  the  so-called  ‘intraseasonal  band.'  Peaks  near 
60  days,  40  days,  and  25  days  are  immediately  evident. 

The  meridional  velocity  spectrum  is  dominated  by 
a  peak  at  ~  360  days,  produced  by  the  annual  vari¬ 
ation  of  the  global  north-south  circulation  (Vincent, 
1993).  A  smaller  semiannual  component  is  also  evi¬ 
dent  here,  which  is  absent  at  lower  heights.  Enhanced 
spectral  power  at  ~  2  days  is  produced  by  a  strong 
quasi-2  day  wave  over  Christmas  Island  (Harris  and 
Vincent,  1993;  Palo  and  Avery,  1996).  However,  in 
contrast  to  the  earlier  study  of  Eckermann  and  Vinc¬ 
ent  (1994)  in  which  no  clear  intraseasonal  peaks  could 
be  identified  in  the  meridional  velocities,  the  longer 
time  series  used  here  reveal  small  but  distinct  peaks 
in  the  meridional-velocity  spectrum  at  ~  60  days  and 
~  35  days. 

Since  a  peak  near  60  days  is  common  to  the  spectra 
of  both  velocity  components  in  Fig.  2,  we  band-pass 
filter  the  amplitude  spectra  using  a  Hanning  (cosine) 
filter  which  is  centered  near  60  days,  as  shown  with  a 
dotted  curve  on  the  spectral  plots  of  Fig.  2.  The  result¬ 
ing  band-passed  oscillations  are  plotted  on  the  bottom 
row  of  Fig.  2. 

Since  60  days  is  approximately  the  third  harmonic 
of  the  strong  MLTSAO  peak  at  —  1 83  days,  we  must 
first  ensure  that  this  quasi-60day  zonal  peak  is  not 
produced  solely  by  higher  harmonics  of  a  non- 
sinusoidal  MLTSAO.  The  20  day  running  averages  on 
the  second  row  of  Fig.  2  are  particularly  useful  in 
assessing  this.  First,  it  is  clear  on  inspecting  these 
smoothed  time  series  that  the  MLTSAO  is  strongly 
nonsinusoidal,  and  so  will  contribute  higher  harmon¬ 
ics.  However,  it  is  also  clear  that  there  are  distinct 
intraseasonal  oscillations  in  the  smoothed  zonal-vel¬ 
ocity  time  series:  see,  for  example,  the  two  westward 
MLTSAO  phases  during  1992. 

Our  new  finding  here  of  similar  (but  smaller)  oscil¬ 
lations  in  the  meridional  velocities  confirms  the  reality 
of  the  quasi-60day  intraseasonal  signals.  The 
smoothed  time  series  of  meridional  wind  in  Fig.  2  also 
shows  clear  intraseasonal  "wiggles,’  this  time  super¬ 
imposed  upon  a  basic  annual  oscillation.  Since  the 
MLTSAO  is,  at  best,  extremely  weak  in  the  meridional 
winds,  spectral  leakage  from  a  nonsinusoidal 


MLTSAO  cannot  explain  the  60  day  meridional  peak 
observed  here.  Note  too  that  the  band-passed  oscil¬ 
lations  of  the  meridional  wind  on  the  bottom-right  of 
Fig.  2  show  bursts  of  amplitude  that  (apart  from  early 
1990)  match  closely  corresponding  bursts  in  the  band- 
passed  zonal  oscillations.  Figure  3  plots  the  20  day 
running  averages  of  the  time  series  at  all  heights  in  the 
range  70-100  km.  These  plots  show  that  the  intra¬ 
seasonal  oscillations  evident  in  the  smoothed  time  ser¬ 
ies  of  Fig.  2  are  coherent  in  both  amplitude  and  phase 
over  a  number  of  heights,  further  reinforcing  the 
reality  of  these  oscillations. 

Thus  clear  intraseasonal  oscillations  of  the  winds  at 
these  heights  are  revealed  by  both  spectral  analysis 
and  inspection  of  the  velocity  time  series.  That  said, 
there  are  two  occasions,  in  early  1993  and  1995,  when 
the  transitions  to/from  strong  westward  MLTSAO 
phases  are  so  rapid  that  some  spectral  power  will  leak 
into  higher  harmonics  of  the  zonal- velocity  spectrum. 
Thus,  some  of  the  quasi-60  day  enhancements  in  the 
band-passed  zonal-velocity  oscillations  observed  at 
these  times  may  be  contributed  to  by  the  MLTSAO. 
Enhancements  observed  at  other  times,  however,  and 
at  all  times  in  the  meridional  time  series,  are  due  to 
distinct  quasi-60  day  oscillations  of  the  winds.  Simi¬ 
larly,  the  peaks  at  ~  40  days  and  20-25  days  in  both 
time  series  cannot  arise  from  a  leakage  of  spectral 
power  from  the  MLTSAO. 

Figure  4  shows  power  spectra  of  the  zonal  and 
meridional  velocity  time  series,  plotted  at  their  heights 
of  observation.  A  corresponding  power  spectral  den¬ 
sity  scale  is  given  on  the  right:  note  that  the  meridional 
spectral  scaling  is  an  order  of  magnitude  smaller  than 
for  the  zonal  spectra.  We  see  that  the  peaks  previously 
identified  at  92  km  in  Fig.  2  occur  over  the  height 
range  80-98  km.  These  spectra  were  also  recomputed 
after  applying  Hanning,  Parzen  (triangular)  and 
Welch  (squared  triangular)  windows  to  the  time  series 
prior  to  Fourier  transformation  to  minimize  any  spec¬ 
tral  leakage.  In  all  cases,  the  major  peaks  identified 
from  Fig.  2  and  evident  in  Fig.  4  occurred  promi¬ 
nently. 

Based  on  these  spectra,  we  choose  three  principal 
pass  bands  (hereafter  labeled  A,  B,  and  C)  which  are 
centered  near  the  major  peaks  in  the  zonal  velocity 
spectra.  These  bands  are  marked  with  rectangles  in 
Fig.  4,  indicating  both  the  width  and  approximate 
height  range  of  the  signals  we  are  trying  to  isolate. 
The  Hanning  filters  we  use  to  isolate  fluctuations 
within  each  band  are  plotted  with  dotted  curves  above 
these  shaded  regions.  The  band-pass  filters  A  and 
C  are  also  applied  to  the  meridional  velocities,  even 
though  there  are  no  distinct  meridional  peaks  within 
band  C.  However,  we  replace  filter  B  when  analyzing 
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Zonal  Winds  at  Christmas  Island  (2°N,  157°W) 


Meridional  Winds  at  Christmas  Island  (2°N,  157°W) 


Fig.  3.  Time-height  plots  of  the  time  series  of  zonal  and  meridional  winds  after  smoothing  with  a  20  day 
running  average.  The  zero  line  of  each  time  series  is  centered  at  the  height  of  observation,  with  the  top 
trace  measured  at  98  km.  Corresponding  velocity  scales  are  shown  on  the  right. 


the  meridional  time  series  with  a  filter  D  which  is 
displaced  to  slightly  smaller  periods,  so  that  it  is  cent¬ 
ered  more  on  the  peak  at  ~35  days  in  the  meridional 
spectra  in  Fig.  4.  Note  that  the  filters  A-C  differ  from 
those  employed  by  Eckermann  and  Vincent  (1994): 
filters  A  and  B  differ  only  slightly,  whereas  filter  C 
isolates  periods  ~  20-30  days  in  this  study,  whereas  in 
Eckermann  and  Vincent  (1994)  it  spanned  periods 
~  15-25  days. 

The  periods  or  period  ranges  of  some  intraseasonal 
phenomena  which  could  conceivably  produce  signals 
in  our  data  are  also  shown  in  Fig.  4.  Recent  studies 
have  shown  evidence  of  the  quasi- 16  day  Rossby  nor¬ 
mal  mode  in  the  mesosphere  at  middle  to  high  lati¬ 
tudes  (Williams  and  Avery,  1992;  Forbes  et  al.,  1995) 


and  in  the  equatorial  ionosphere  (Forbes  and 
Leveroni,  1992;  Parish  et  al.,  1995).  Modeling  by 
Forbes  et  al.  (1995)  suggests  that  the  penetration  of 
this  disturbance  into  the  equatorial  mesosphere  and 
thermosphere  depends  sensitively  on  background 
wind  profiles  and  the  effects  of  gravity-wave  drag. 
There  is  no  clear  evidence  of  a  16  day  signal  in  Fig.  4 
(although  there  is  a  smattering  of  activity  throughout 
the  12-20  day  band  within  which  this  mode  sometimes 
resides),  nor  is  their  much  evidence  of  16  day  Rossby 
modes  in  the  equatorial  troposphere  and  stratosphere 
(Hirooka  and  Hirota,  1989).  Clearly  a  quasi-16day 
wave  cannot  account  for  the  major  peaks  in  Fig.  4. 
Similarly,  lunar  semidiurnal  tides,  while  evident  in 
mesospheric  wind  data  from  extratropical  sites  (e.g., 
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Stening  et  al.,  1994),  do  not  appear  prominently  in 
these  data  either:  indeed,  Stening  et  al.  (1997)  found 
amplitudes  of  no  more  than  l-2ms~‘  in  these  data. 
Neither  is  there  much  evidence  for  a  response  to  the 
27  day  variation  of  incoming  ultraviolet  radiation  due 
to  solar  rotation,  which  is  known  to  influence  meso¬ 
spheric  ozone  and  temperatures,  particularly  near  the 
equator  (e.g.,  Summers  et  al.,  1990;  Hood  et  al.,  1991; 
Brasseur,  1993).  Indeed,  the  only  phenomena  that 
appear  to  relate  to  the  peaks  we  observe  are  the  30- 
60  day  and  20-25  day  oscillations  of  the  equatorial 
troposphere  (e.g.,  Hartmann  et  al.,  1992;  Madden  and 
Julian,  1994),  as  argued  by  Eckermann  and  Vincent 
(1994). 


TIME-HEIGHT  INTRASEASONAL  VARIABILITY 

Figure  5  profiles  the  zonal  velocity  activity  within 
various  period  bands.  The  top  plot  shows  the  mean 
winds,  obtained  by  retaining  periods  >  100  days.  The 
semiannual  oscillation  is  clearly  evident,  showing,  for 
example,  that  the  first  westward  phase  of  the 
MLTSAO  is  usually  stronger  (e.g.,  Garcia  et  al., 
1997).  Considerable  interannual  variability  is  also  evi¬ 
dent:  for  example,  there  are  atypically  strong  first 
westward  phases  in  early  1993  and  1995,  whereas  both 
westward  phases  are  weaker  and  roughly  equal  during 
1992  and  1997.  Note  also  the  much  smaller  semi¬ 
annual  oscillation  below  ~75km.  It  is  well-known 
that  the  amplitude  of  the  semiannual  oscillation  is 
small,  between  about  60-70  km,  effectively  separating 
the  stratospheric  semiannual  oscillation  (SSAO)  and 
the  MLTSAO  (e.g.,  Hirota,  1978;  Garcia  et  al.,  1997). 
Band-passed  fluctuations  within  bands  A-C  are  dis¬ 
played  in  the  remaining  three  plots  in  Fig.  5,  and 
again,  the  region  below  ~75km  is  relatively  quiesc¬ 
ent,  whereas  at  80-98  km  there  are  transient  bursts  of 
activity  within  all  three  bands. 

Comparisons  of  Christmas  Island  MF  radar  winds 
with  coincident  wind  measurements  from  UARS  have 
suggested  significant  underestimations  of  the  radar 
winds  at  lower  heights,  but  much  better  agreement 
higher  up  (see  Figs  1 1  and  12  of  Burrage  et  al.  (1996)). 
At  these  lower  heights,  the  radar’s  data  rates  decrease, 
diurnal  data  coverage  is  incomplete,  and  clutter  from 
ocean-wave  echoes  becomes  a  problem  (Vincent  and 
Lesicar,  1991),  all  of  which  could  contribute  to  these 
apparent  underestimations  of  wind  speeds.  Thus,  the 
absence  of  intraseasonal  peaks  and  the  small 
MLTSAO  below  ~  75  km  may  also  be  due  in  part  to 
degraded  instrument  sensitivity  at  these  lower  heights. 

The  bursts  of  activity  within  bands  A  and  B  gen¬ 
erally  occur  at  similar  times.  As  noted  by  Eckermann 


and  Vincent  (1994),  the  occurrence  of  activity  in  these 
bands  is  well  correlated  with  the  strength  of  the  west¬ 
ward  MLTSAO  zonal  winds.  Since  the  first  westward 
phase  is  usually  stronger  than  the  second,  there  tends 
to  be  more  activity  in  bands  A  and  B  at  the  start  of 
the  year.  Activity  in  band  C,  on  the  other  hand,  tends 
to  appear  in  the  middle  to  latter  half  of  the  year. 
Active  regions  within  all  three  bands  descend  phase 
with  time,  in  agreement  with  the  initial  study  of  bands 
A  and  B  by  Eckermann  and  Vincent  (1994). 

To  study  possible  seasonal  variations  in  the  intra¬ 
seasonal  activity,  we  apply  phase-time  methods  to  the 
band-passed  time  series  (e.g.,  Huang  et  al.,  1992).  The 
basic  method  is  depicted  in  Fig.  6.  Its  advantage  here 
is  that  it  provides  time  series  of  the  peak  amplitude 
A(t)  and  phase  4>(t)  of  the  intraseasonal  signals  at 
the  daily  resolution  of  the  original  time  series,  unlike 
variance  calculations  which  degrade  the  resolution  to 
intraseasonal  time  scales  (unless  some  sliding  average 
is  employed)  and  make  seasonal  trends  harder  to 
identify.  Figure  7  displays  seasonal  variations  of  the 
peak  amplitudes  of  the  band-passed  time  series  in  Fig. 
5.  Despite  the  interannual  variability  evident  within 
all  three  bands,  we  note  weak  but  distinct  annual 
variability  in  the  mean  seasonal  curves.  Band  A  shows 
greatest  activity  around  December-April,  and  a  simi¬ 
lar  tendency  is  weakly  evident  in  band  B.  Conversely, 
band  C  shows  greater  activity  on  average  during 
August-December. 

Phase-time  methods  can  also  provide  time  series  of 
frequency  co(t)  (see  Fig.  6),  and  hence  the  periods 
of  various  intraseasonal  bursts  of  activity.  Figure  8 
profiles  height  variations  of  the  period  and  peak 
amplitude  of  four  bursts  of  band-A  activity  that 
occurred  early  in  1990,  1991,  1993,  and  1995  (see  Fig. 
5).  Activity  typically  peaks  at  ~  86  km,  but  attenuates 
rapidly  below  80  km  as  previously  noted.  Periods  in 
these  active  regions  are  in  the  range  55-65  days,  vary¬ 
ing  slightly  from  event  to  event,  although  the  results 
also  depend  somewhat  on  the  filter  functions  used. 
Similar  calculations  for  four  burst  regions  in  band  B 
(during  early  1990,  early  1992,  late  1992,  and  early 
1994)  showed  similar  characteristics,  with  peak  ampli¬ 
tudes  ~5-8ms~'  and  periods  ~  40  days.  Four  burst 
regions  in  the  band-C  time  series  (in  mid- 1990,  mid- 
1991,  mid- 1992,  and  late  1993)  showed  amplitudes 
peaking  at  ~5-8ms~'  near  ~ 85-90  km,  attenuated 
amplitudes  below  80  km,  and  periods  typically  ~  22- 
24  days. 

Figure  9  profiles  the  activity  in  the  meridional  vel¬ 
ocities.  The  variations  in  the  band-A  activity  are  simi¬ 
lar  to  those  in  Fig.  5,  and  these  active  regions  also 
exhibit  phase  descent  with  time.  The  peak  of  the 
activity  tends  to  occur  typically  ~90km,  slightly 
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Fig.  4.  Power  spectral  densities  of  zonal  and  meridional  velocities  within  the  intraseasonal  band,  with  the 
zero  line  plotted  at  the  height  of  observation.  Spectral  units  are  given  on  the  right.  Bands  A-C  are  marked 
on  these  plots,  and  the  Hanning  filters  used  to  isolate  these  peaks  are  plotted  with  dotted  curves  above  the 
marked  areas.  Periods  and  period  ranges  of  other  intraseasonal  geophysical  phenomena  are  depicted  on 

the  plot. 
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Zonal  Wind  (Periods  >108  days)  over  Christmas  Island  (2°N,  157°W) 
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Fig.  5.  Time-height  plots  of  the  band-passed  time  series  of  zonal  winds  at  various  heights  are  shown  in  the 
bottom  three  plots  for  filters  A-C.  An  amplitude  scale  for  these  time  series  is  shown  on  the  right.  The 
activity  in  these  oscillations  can  be  compared  with  activity  in  the  top  plot,  which  shows  contours  of  zonal 
winds  in  which  periods  >  100  days  have  been  retained,  which  we  take  to  represent  mean-flow  variations. 
Contour  labels  are  in  ms-1,  and  westward  (eastward)  values  are  negative  (positive)  and  have  darker 
(lighter)  shading  and  solid  white  (black)  contours.  The  contour  interval  is  10ms-1.  Zero  wind  lines  are 

dotted. 
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Fig.  6.  Diagram  showing  the  basis  of  the  phase-time  method. 
The  original  time  series  U(t)  is  combined  with  its  Hilbert 
transform,  which  we  approximate  here  by  applying  90  phase 
shifts  to  each  harmonic  of  the  Fourier  transform  of  U(t), 
then  retransforming  to  yield  U+W(t).  The  two  time  series  are 
combined  into  a  complex  vector  time  series  (solid  vector) 
which  provides  both  the  peak  amplitude  A(t )  (the  length  of 
Ihe  solid  vector)  and  the  phase  <j>(t).  A  time  series  of  frequency 
follows  from  w(t)  =  d<j>(t)!dt. 


higher  than  for  the  zonal  activity.  The  major  differ¬ 
ence,  however,  is  that  the  meridional  oscillations  in 
band  A  are  out  of  phase  with  the  zonal  band-A  oscil¬ 
lations.  This  is  shown  in  Fig.  10a,  which  plots  time 
hodographs  of  the  total  (zonal  and  meridional)  band- 
A  velocities  at  90  km.  An  anticorrelation  between  the 
band-A  velocity  components  is  apparent,  with  cor¬ 
relation  coefficients  of  around  — 0.75  at  this  height 
over  the  full  1944  days  of  data  (Fig.  10b). 

CORRELATIONS  WITH  GRAVITY  WAVE  AND  TIDAL 
ACTIVITY 

Figure  1 1  plots  time  series  of  the  zonal  and  mer¬ 
idional  winds  from  Fig.  3,  along  with  the  peak  zonal 
amplitude  of  the  diurnal  tide,  and  zonal  and  mer¬ 
idional  velocity  variances  in  the  12min^lh  band,  the 
latter  indicative  of  gravity-wave  activity.  At  times 
when  intraseasonal  activity  was  strong  (e.g.,  early  in 
1990,  1991,  and  1993),  we  see  intraseasonal  oscil¬ 
lations  of  zonal  gravity-wave  variances  and  diurnal 
tidal  amplitudes  above  ~  80-85  km  which  are  coher¬ 
ent  in  phase  with  height.  The  meridional  velocity  vari¬ 
ances  tend  to  show  strongest  intraseasonal  oscillations 
above  about  90  km. 

To  investigate  the  gravity- wave  variability  further, 
we  compute  power  spectra  of  the  unsmoothed  time 
series  of  zonal  and  meridional  velocity  variances,  to 
delineate  any  periodicities  in  the  wave  activity  (e.g., 
Isler  and  Fritts,  1996;  Connor  and  Avery,  1996).  The 
spectra  at  86-90  km  are  shown  in  Fig.  12a.  A  strong 
peak  at  ~  60  days  is  evident  at  86-90  km  in  the  zonal 
variances,  which  clearly  is  very  similar  to  the  quasi- 


60  day  band-A  peak  in  the  zonal  wind  spectra  from 
Fig.  4.  The  intraseasonal  variations  at  86-90  km  are 
greater  for  the  zonal  variance  than  for  the  meridional 
variance.  There  are  suggestions  of  additional  period¬ 
icities  in  both  velocity  variances  at  ~  40  days  (band 
B),  and  ~  25  days  (band  C),  although  these  peaks  are 
quite  weak  and  may  not  be  significant. 

The  variance  spectra  above  86-90  km  change  mark¬ 
edly.  In  Fig.  12b  we  note  that  intraseasonal  variations 
in  zonal  variance  at  92-98  km  are  of  similar  magnitude 
to  those  at  86-90  km.  The  peak  at  60  days,  however, 
has  been  strongly  dissipated  (along  with  the  weak 
40  day  and  25  day  peaks),  consistent  with  a  transfer 
of  these  periodicities  in  gravity-wave  activity  into  the 
mean  flow  through  the  86-98  km  height  range,  as 
hypothesized  by  Eckermann  and  Vincent  (1994).  The 
zonal  variance  spectrum  at  92-98  km  now  peaks 
nearer  45  days.  Intraseasonal  variations  of  meridional 
activity  at  92-98  km  are  large  and  exceed  the  zonal 
variations.  A  strong  peak  in  the  meridional  variances 
now  occurs  at  ~  60  days,  with  subsidiary  peaks  at 
~  90  days  and  ~  45  days. 

Figure  13  shows  spectra  of  time  series  of  the  peak 
zonal  and  meridional  amplitudes  of  the  diurnal  tide. 
Again,  a  quasi-60day  peak  is  evident.  However,  the 
zonal  tidal  amplitudes  also  show  large  distinct  peaks 
centered  at  ~  40  days  and  ~  24  days,  which  are  very 
similar  to  the  peaks  within  bands  B  and  C  of  the  zonal- 
wind  spectra  in  Fig.  4  and  far  exceed  the  weak  peaks 
in  the  gravity-wave  variances  in  Fig.  12a. 

The  intraseasonal  tidal  activity  and  the  attenuation 
(intensification)  of  intraseasonal  variations  of  zonal 
(meridional)  gravity-wave  activity  with  height  are  all 
well  illustrated  during  the  strong  MLTSAO  phase  in 
early  1993.  The  winds,  zonal  tidal  amplitudes,  and 
gravity-wave  activity  for  the  first  six  months  of  1993 
are  replotted  in  Fig.  14.  We  see  strong  intraseasonal 
variations  in  zonal  gravity-wave  variance  which  atten¬ 
uate  with  height,  consistent  with  the  perceived  role  of 
this  gravity-wave  activity  in  driving  both  the  strong 
westward  MLTSAO  winds  and  intraseasonal  wind 
variability.  We  also  note  a  very  strong  diurnal  tide 
with  an  intraseasonal  activity  variation  which  atten¬ 
uates  with  height,  which  may  also  drive  these  strong 
westward  winds.  The  strong  colocated  gravity- wave 
and  tidal  activities  also  produce  strong  gravity-wave 
tidal  interactions  in  these  data  (see  also  Nakamura  et 
al.,  1997).  The  intensification  of  meridional  variance 
at  upper  heights  is  also  very  apparent,  as  is  its  intra¬ 
seasonal  variation.  We  tentatively  associate  this  with 
the  filtering  of  more  meridionally  directed  waves  by 
the  intraseasonal  flow  variations  that  appear  at  ~  80- 
90  km.  Note  the  quasi-bimonthly  variations  in  all  five 
plots. 
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Seasonal  Variation  of  Intraseasonal  Amplitudes 
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Fig.  7.  Seasonal  variations  of  peak  amplitudes  of  the  band-passed  zonal  wind  time  series  profiled  in  Fig. 
5.  Various  broken  curves  show  the  trends  for  individual  years,  while  the  thick  curve  shows  the  mean 

seasonal  variation  from  1990-1994. 


discussion  ~  60  days  (band  A),  35-40  days  (bands  B  and  D),  and 

20-25  days  (band  C).  Similar  periodicities  have  also 
been  detected  in  diurnal  tidal  amplitudes  and  gravity- 
wave  variances  at  these  heights.  We  have  also  pre- 


We  have  identified  three  distinct  bands  of  intra¬ 
seasonal  activity  in  these  MLT  wind  data,  centered  at 
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Band  A:  Day  Numbers  395-500. 
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Fig.  8.  Altitude  profiles  of  period  (left  column)  and  peak  zonal  velocity  amplitude  (right  column)  for  four 
separate  bursts  of  activity  within  band  A.  The  day  number  ranges  are  shown  on  the  plot,  where  day  1  is 
January  1.  1990.  Error  bars  are  standard  errors  of  the  mean  derived  from  the  time  series  of  amplitude  and 

phase  provided  by  phase-time  methods. 
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BAND  C:  Meridional  Wind  Oscillations 
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Fig.  9.  Analogous  presentation  to  Fig.  5,  but  this  time  for  meridional  winds.  The  amplitude  scale  for  the 
band-passed  fluctuations  is  shown  on  the  right  of  the  bottom  three  plots.  The  contour  interval  in  the  top 
plot  is  5  ms-1,  and  southward  values  are  negative  and  have  darker  shading  and  white  solid  contours.  Note 
that  filter  B  is  replaced  by  filter  D  in  this  case. 
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(a)  Band  A  Velocity  Hodograph  at  z  =  90  km. 
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Fig.  10.  (a)  Hodograph  of  the  band-A  velocity  oscillations  at  90  km.  (b)  Height  variation  of  the  correlation 
coefficient  between  the  zonal  and  meridional  velocity  time  series  within  band  A. 


sented  preliminary  evidence  (summarized  in  Fig.  4) 
that  these  intraseasonal  winds  may  be  associated  with 
the  30-60  day  and  20-25  day  oscillations  encountered 
in  the  equatorial  troposphere,  as  postulated  by  Eck- 
ermann  and  Vincent  (1994). 

We  investigate  this  in  more  depth  here.  Since  the 
data  from  Christmas  Island  comprise  our  only  explicit 
information  on  intraseasonal  oscillations  of  the  equa¬ 
torial  MLT  region  to  date,  firm  conclusions  on  the 
nature  and  origin  of  this  activity  are  difficult  to  draw. 
Our  discussion  is  developed  with  this  limitation  in 
mind,  and  is  structured  as  follows.  First  we  briefly 
describe  the  salient  characteristics  of  observed  tropical 


tropospheric  intraseasonal  variability  (TTIV).  Next, 
based  on  some  recent  observations  in  and  around 
India,  we  briefly  reconsider  whether  planetary  waves 
associated  with  TTIV  could  propagate  directly  to 
mesospheric  heights.  Finally,  on  the  weight  of  current 
evidence,  we  propose  an  explanation  for  the  intra¬ 
seasonal  oscillations  we  observe  based  on  gravity- 
wave  and  diurnal-tidal  driving  of  the  MLT  flow,  and 
go  on  to  investigate  some  possible  consequences. 

Tropical  tropospheric  intraseasonal  variability  (TTIV) 

The  30-60  day  tropospheric  oscillation  is  an  inten¬ 
sively  studied  phenomenon,  and  is  known  as  the  Mad- 
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Fig.  1 1.  The  top  two  plots  reproduce  time  series  from  1990-1994  of  zonal  and  meridional  winds  from  Fig. 
3.  The  middle  plot  shows  peak  zonal  amplitudes  of  the  diurnal  tide.  The  bottom  two  plots  show  gravity- 
wave  variances  in  zonal  and  meridional  velocities  within  the  1 2  rnin-4  h  period  range.  The  gravity- wave 
time  series  were  smoothed  with  a  10  day  running  average,  and  the  tidal  data  with  a  6  day  running  average. 
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Power  Spectra  of  Velocity  Variance  (12  min-4  hours) 


(a)  86-90  km 
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Fig.  12.  Power  spectra  of  time  series  of  zonal  and  meridional  velocity  variance  in  the  12min-4h  band: 
(a)  86-90  km;  (b)  92-98  km.  These  spectra  were  smoothed  with  a  3-point  running  average. 


den-Julian  oscillation  (MJO),  after  its  discoverers  (see 
the  review  of  Madden  and  Julian  (1994)).  It  occurs 
when  a  large  region  of  convective  ascent  moves  east¬ 
ward  from  the  Indian  Ocean  into  the  western  Pacific 
(Madden  and  Julian,  1972),  producing  30-60 day 
periodicities  in  regional  convective  activity  as  it  passes 
(see,  e.g.,  Hartmann  and  Gross  1988;  Salby  and 
Hendon,  1994;  Dunkerton  and  Crum,  1995).  This  dis¬ 
turbs  the  climatological  Walker  circulation,  pro¬ 
ducing  a  planetary-scale  response  of  similar  period  in 
the  equatorial  velocity  field.  This  velocity  response 
can  be  understood  to  first  order  as  a  coupled  Kelvin- 
Rossby  mode  near  the  forcing  region,  which,  as  one 
shifts  further  eastward,  is  replaced  by  a  wavenumber- 
1  Kelvin  wave  radiating  away  from  this  forcing  zone 


(e.g.,  Parker,  1973;  Hendon  and  Salby,  1994,  1996; 
Salby  et  al.,  1994). 

A  strong  20-25  day  tropical  oscillation  has  also 
been  detected  over  the  western  Pacific  (Ghil  and  Mo, 
1991;  Hartmann  et  al.,  1992;  Kiladis  et  al.,  1994)  and 
elsewhere  (e.g.,  Anyamba  and  Weare,  1995).  While 
this  particular  oscillation  seems  to  be  distinct  from  the 
MJO,  it  should  also  be  noted  that  the  MJO's  period 
range  can  vary  considerably,  and  occasionally 
approaches  25  days  (e.g.,  Hayashi  and  Golder,  1993; 
Madden  and  Julian,  1994). 

Can  TTIV propagate  into  the  middle  atmosphere? 

Eckermann  and  Vincent  (1994)  argued  on  theor¬ 
etical  grounds  that  the  mesospheric  intraseasonal 
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Power  Spectra  of  Diurnal  Tidal  Amplitudes 
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Fig.  13.  Power  spectra  of  time  series  of  the  peak  zonal  (solid 
curve)  and  meridional  (dotted  curve)  velocity  amplitudes  of 
the  diurnal  tide:  (a)  86-92  km;  (b)  92-98  km.  These  spectra 
were  smoothed  with  a  3-point  running  average. 


activity  observed  at  Christmas  Island  was  unlikely  to 
result  from  direct  upward  propagation  of  the  plan¬ 
etary  waves  associated  with  the  MJO.  However,  some 
recent  observations  in  the  Indian  sector  lead  us  to 
reconsider  this  possibility  here. 

It  is  well  known  that  MJO  activity  attenuates  rap¬ 
idly  above  the  tropopause  (Madden  and  Julian,  1971), 
and  studies  of  rocket  and  radiosonde  data  front  trop¬ 
ical  Indian  stations  show  this  quite  clearly  (Nagpal 
and  Raghavarao,  1991;  Nagpal  et  al.,  1994;  Kumar 
and  Jain,  1994).  However,  Nagpal  and  Raghavarao 
(1991)  showed  that  intraseasonal  power  in  their  rocket 
wind  data  reintensified  somewhat  in  the  upper  strato¬ 
sphere.  This  was  investigated  further  by  Nagpal  et  al. 
(1994)  and  Kumar  and  Jain  (1994)  using  rocket  data 
from  three  Indian  stations  spanning  8.5°-21.5°N 
(shown  in  Fig.  1).  They  noted  increased  intraseasonal 
activity  in  both  zonal  and  meridional  winds  at 
~  50  km,  most  notably  at  the  two  higher-latitude  sites. 
Since  there  were  nearly  equal  amounts  of  intra¬ 


seasonal  activity  in  both  wind  components,  they  con¬ 
cluded  that  this  was  inconsistent  with  Kelvin  waves, 
and  associated  the  activity  with  Rossby  waves.  Nagpal 
and  Raghavarao  (1991)  argued  that  these  Rossby 
waves  propagated  into  the  tropical  upper  stratosphere 
from  mid-latitudes  (see  also  Nagpal  et  al.  (1994)), 
whereas  Kumar  and  Jain  (1994)  suggested  that  the 
waves  ‘leaked’  into  the  upper  stratosphere  from  the 
tropical  troposphere  via  vertical  propagation.  The  lat¬ 
ter  is  possible  since  the  MJO  contains  a  significant 
Rossby-wave  component,  particularly  in  forcing 
zones  near  India  (e.g.,  Hendon  and  Salby,  1994). 

Ziemke  and  Stanford  (1991)  studied  Rossby  waves 
of  intraseasonal  period  (1-2  months)  in  eight  years 
of  daily  global  geopotential  height  data  from  British 
Meteorological  Office  analyses.  They  focused  on 
strong  wave  activity  in  a  longitude  zone  near  the  east 
coast  of  India.  This  activity,  initially  concentrated  in 
the  equatorial  troposphere,  propagated  into  the  extra- 
tropical  upper  troposphere  and  stratosphere  to  form 
a  regional  dipole  pattern  about  the  equator.  Similar 
patterns  are  also  observed  at  periods  ~  25  days 
(Anyamba  and  Weare,  1995).  Thereafter,  the  waves 
continued  propagating  upwards,  refracting  back  into 
the  tropical  stratopause  region  (see  Figs  5  and  6  of 
Ziemke  and  Stanford  (1991)).  Similar  Rossby-wave 
propagation  paths  are  suggested  in  the  model  simu¬ 
lations  of  Salby  et  al.  (1994)  (e.g.,  their  Fig.  18).  This 
sort  of  wave  propagation  path  could  clearly  explain 
the  reemergence  of  intraseasonal  activity  observed  in 
the  Indian  rocket  data  near  the  stratopause:  if  so, 
then  the  explanations  of  it  offered  by  Nagpal  and 
Raghavarao  (1991)  and  Kumar  and  Jain  (1994)  are 
essentially  both  correct. 

These  rocket  and  satellite  data,  however,  cut  out  at 
~60km,  and  so  it  is  unclear  whether  these  waves 
reach  the  mesosphere  and  lower  thermosphere.  GCM 
simulations  by  Hayashi  and  Golder  (1993)  show  pen¬ 
etration  of  equatorial  Rossby-wave  activity  at  intra¬ 
seasonal  periods  to  near  the  top  of  their  model 
(~75km),  which  is  still  below  the  region  of  activity 
observed  in  our  data  (~  80-98  km).  However,  the 
presence  of  Rossby  waves  at  the  tropical  stratopause 
is  only  feasible  during  the  eastward  phases  of  the 
SSAO.  Since  the  antiphased  MLTSAO  will  be  west¬ 
ward  at  these  times,  equatorial  Rossby  waves  will 
be  absorbed  at  critical  levels  somewhere  between  the 
SSAO  and  MLTSAO.  Thus  direct  Rossby-wave 
propagation  into  the  ~  80-98  km  region  of  the  equa¬ 
torial  atmosphere  seems  unlikely  at  these  times.  How¬ 
ever,  any  intraseasonal  Rossby-wave  oscillations  at 
the  stratopause  could  affect  the  transmissivity  of  grav¬ 
ity  waves  through  the  tropical  stratopause,  thereby 
forcing  a  mesospheric  disturbance  of  similar  period 
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Fig.  14.  As  for  Fig.  11.  but  focusing  on  the  first  six  months  of  1993,  with  the  data  plotted  in  contour  form. 
Contour  labels  are  in  ms-1  for  the  mean  winds  and  zonal  diurnal  tidal  amplitudes,  and  in  m2s-’  for  the 
velocity  variances.  Active  regions  in  Ihe  data  are  highlighted  with  heavier  contours. 
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through  modulated  gravity-wave  breaking.  Very  simi¬ 
lar  processes  have  been  proposed  in  the  mid-latitude 
Northern  Hemisphere  during  winter  (Holton,  1984; 
McLandress  and  McFarlane,  1993;  Smith,  1996). 

Intraseasonal  activity  in  the  mid-Pacific  MLT  region 

Our  mid-Pacific  observations  at  MLT  heights  have 
shown  that  intraseasonal  activity  within  bands  A  and 
B  occurs  preferentially  during  westward  phases  of  the 
MLTSAO.  As  mentioned,  any  equatorial  Rossby 
waves  have  probably  been  removed  at  critical  levels 
lower  down,  but  since  the  SSAO  is  eastward  at  these 
times,  intraseasonal  Rossby  waves  could  be  present  at 
the  stratopause  to  modulate  gravity-wave  trans¬ 
missivity  and  force  a  related  mesospheric  disturbance 
(Holton,  1984).  Kiladis  and  Wheeler  (1995)  have 
observed  Rossby  waves  with  periods  ~  6-30  days  in 
the  troposphere  near  Christmas  Island,  although 
whether  these  waves  could  propagate  into  the  tropical 
upper  stratosphere  is  unclear.  GCM  simulations  by 
Hayashi  and  Golder  (1993)  produced  significant  intra¬ 
seasonal  Rossby-wave  activity  throughout  the  equa¬ 
torial  stratopause. 

However,  there  are  several  reasons  for  supposing 
that  much  of  the  activity  at  80-98  km  over  Christmas 
Island  is  not  produced  either  directly  or  indirectly 
by  middle-atmosphere  Rossby  waves.  First,  Rossby- 
wave  oscillations  or  Rossby-wave-modulated  gravity- 
wave  breaking  might  be  expected  to  accelerate  zonal 
and  meridional  winds  fairly  equally  (see  e.g.,  Kumar 
and  Jain,  1994),  whereas  the  bulk  of  the  intraseasonal 
activity  we  observe  here  is  zonal.  The  activity  over 
Christmas  Island  is  more  suggestive  of  either  a  zonal- 
wind  vacillation  or  of  a  Kelvin  wave  (Eckermann  and 
Vincent,  1994).  Indeed,  the  data  of  Ziemke  and  Stan¬ 
ford  (1991)  suggest  that  the  upper-stratospheric 
Rossby  waves  they  studied  were  confined  to  a  fairly 
narrow  longitude  zone  near  India,  and  so  may  be  less 
relevant  to  our  observations  in  the  central  Pacific. 

This  perceived  longitudinal  difference  in  the  charac¬ 
ter  of  intraseasonal  activity  at  the  Indian  stratopause 
and  in  the  mid-Pacific  MLT  region  has  interesting 
similarities  with  the  MJO.  In  the  forcing  region  in 
and  around  Indonesia  and  India,  MJO  activity  is 
characterized  by  a  combined  Kelvin-Rossby  wave 
response,  whereas  at  longitudes  well  away  from  the 
forcing  region  the  activity  is  dominated  by  a  pure 
(zonal)  Kelvin  wave  response  (Hendon  and  Salby, 
1994,  1996).  Tropospheric  radiosonde  data  from  Kan- 
ton  (formerly  Canton  Island),  which  is  very  near  Chri¬ 
stmas  Island  (see  Fig.  1),  bear  this  out.  They  reveal  a 
strong  peak  in  the  zonal  wind  at  ~  40-50  days 
(Madden  and  Julian,  1971;  Hartmann  and  Gross, 


1988)  produced  by  a  large-amplitude  Kelvin  wave 
which  has  propagated  eastward  to  Kanton  from 
its  forcing  region  to  the  west  (Hendon  and  Salby, 
1994). 

The  MJO-related  Kelvin  waves  over  Kanton  atten¬ 
uate  rapidly  above  ~80hP  (Madden  and  Julian, 
1971),  supporting  theoretical  assessments  that  they 
cannot  propagate  vertically  to  MLT  heights  (Eck¬ 
ermann  and  Vincent,  1994).  Since  they  are  equa- 
torially  ducted,  an  extratropical  propagation  path  into 
the  upper  atmosphere,  like  the  one  chronicled  for 
Rossby  waves  by  Ziemke  and  Stanford  (1991),  is  also 
infeasible.  How,  then,  might  similar  signals  reach  the 
mid-Pacific  mesosphere? 

Gravity  waves,  tides,  and  TTIV-influenced  convective 
activity 

Eckermann  and  Vincent  (1994)  suggested  that 
MJO-modulated  convection  could  produce  an  intra¬ 
seasonal  variation  in  the  intensity  of  convectively  gen¬ 
erated  gravity  waves  entering  the  middle  atmosphere. 
Models  have  shown  how  the  dissipation  of  con¬ 
vectively  generated  waves  can  greatly  influence  the 
circulation  of  the  middle  atmosphere  (e.g.,  Pfister  et 
al.,  1993a;  Alexander,  1996).  Selected  aircraft  missions 
during  the  Panama  80  and  STEP  87  campaigns  (in 
regions  shown  in  Fig.  1)  revealed  strong  stratospheric 
gravity-wave  activity  emanating  from  tropical  tro¬ 
pospheric  convective  systems  (Pfister  et  al.,  1993a,  b; 
Alexander  and  Pfister,  1995). 

Here  we  have  noted  correlations  between  intra¬ 
seasonal  variations  in  MLT  winds  and  similar  period¬ 
icities  in  gravity-wave  activities  at  80-98  km.  The 
relations  at  86-90  km  are  superficially  consistent  with 
a  60  day  MJO-controlled  variation  in  the  intensity  of 
convectively  generated  gravity  waves  produced  in  the 
troposphere  (Fig.  12).  However,  similar  cycles  in  wave 
variance  could  conceivably  have  resulted  from  gravity 
waves  interacting  with  a  preexisting  quasi-60  day  wind 
oscillation  at  these  heights  (e.g.,  Vincent  and  Ecker¬ 
mann,  1990;  Isler  and  Fritts,  1996). 

Some  evidence  for  the  hypothesized  MJO  modu¬ 
lation  of  gravity-wave  generation  by  convection  in  the 
tropical  Pacific  was  provided  recently  by  Karoly  et  al. 
(1996).  They  analyzed  several  months  of  high-res¬ 
olution  radiosonde  data  from  Santa  Cruz  (see  Fig.  1) 
in  the  tropical  south  Pacific.  It  is  well-known  that 
convective  activity  in  and  around  this  region  is  influ¬ 
enced  by  the  MJO  (see  e.g.,  Hartmann  and  Gross, 
1988)  as  eastward -moving  intraseasonal  convective 
zones  drift  southward  towards  the  South  Pacific  Con¬ 
vergence  Zone  (e.g.,  Kiladis  and  Wheeler,  1995).  The 
observations  of  Karoly  et  al.  (1996)  were  taken  during 
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the  4  month  TOGA  COARE  observational  campaign, 
starting  in  November,  1992  when  the  MJO  was  active 
and  dominated  regional  convective  activity  (Gutzler  et 
al.,  1994;  Sheu  and  Liu,  1995).  They  noted  statistically 
significant  correlations  between  stratospheric  gravity- 
wave  activity  in  their  data  and  outgoing  longwave 
radiation  in  the  region,  the  latter  a  standard  proxy  for 
convective  activity.  Thus,  during  this  period,  there 
appear  to  be  clear  MJO-related  influences  on  the 
intensity  of  gravity  waves  entering  the  middle  atmo¬ 
sphere. 

We  have  seen  even  clearer  intraseasonal  period¬ 
icities  in  the  amplitude  of  the  diurnal  tide  (Fig.  13). 
Convective  activity  in  the  equatorial  troposphere  has 
a  strong  diurnal  variation  (e.g.,  Hendon  and  Wood- 
berry,  1993)  which,  through  latent  heat  release,  can 
provide  an  important  thermal  forcing  of  tidal  struc¬ 
tures  (Hamilton,  1981).  Modeling  has  shown  that 
diurnally  varying  convective  activity  forces  a  rich 
spectrum  of  nonmigrating  diurnal  tidal  modes  in  the 
atmosphere  (Lieberman  and  Leovy,  1995;  Williams 
and  Avery,  1996),  and  the  largest-scale  modes  can 
propagate  into  the  equatorial  MLT  region  with  sig¬ 
nificant  amplitudes  (Lieberman,  1991;  Lieberman  and 
Leovy,  1995).  Given  the  intraseasonal  variations  in 
equatorial  convective  activity  that  also  occur  (e.g., 
Hartmann  and  Gross,  1988),  then  it  seems  likely  that 
combined  diurnal  and  intraseasonal  variations  in  con¬ 
vection  will  yield  intraseasonal  variations  in  the  ampli¬ 
tude  of  convectively  forced  diurnal  tides  in  the  MLT 
region,  as  we  observed  in  Fig.  13. 

Eckermann  and  Vincent  (1994)  argued  that  MJO- 
modulated  gravity-wave  activity  might  propagate  into 
the  mesosphere  and  dissipate,  forcing  similar  MJO- 
related  periodicities  in  mesospheric  winds.  The  same 
argument  also  holds  for  diurnal  tides,  which  dissipate 
at  equatorial  MLT  heights  and  can  provide  sub¬ 
stantial  mean-flow  forcing  (e.g.,  Hamilton,  1995). 
Indeed,  comparison  of  Figs  12  and  13  suggests  that 
diurnal  tides  may  be  more  important  than  gravity 
waves  in  forcing  the  intraseasonal  zonal-wind  activity 
within  bands  B  and  C  (Fig.  4).  The  intraseasonal  vari¬ 
ations  we  have  observed  in  gravity-wave  activity,  tidal 
activity,  and  winds  in  these  MLT  data,  coupled  with 
the  lower-atmosphere  observations  of  Karoly  et  al. 
(1996),  provide  some  preliminary  observational  sup¬ 
port  for  gravity/tidal-wave  forcing  of  intraseasonal 
MLT  winds  near  the  equator.  We  now  develop  these 
ideas  in  a  little  more  theoretical  depth,  and  then  apply 
them  to  investigate  how  TTIV-controlled  gravity- 
wave  and  tidal  activity  could  produce  some  of  the 
observed  intraseasonal  periodicities  we  see  in  the 
MLT  winds. 


Intraseasonal  gravity-wave  and  tidal  forcing  of  winds  in 
the  equatorial  MLT  region 

Near  the  equator,  the  zonally  averaged  momentum 
budget  can  be  described  quite  accurately  by  the  sim¬ 
plified  equation  (e.g.,  Bacmeister  et  a!.,  1995) 

DU  _ 

,  A, (1) 

where  D/Dt  =  d/dt  +  v*dldy  +  w*8/8z,  ( v *,  w*)  are 
the  velocities  of  the  residual  meridional  circulation,  U 
is  zonal  velocity  and  Xmech  is  a  mechanical  forcing  term 
which,  at  MLT  heights,  is  dominated  by  gravity  waves 
and  tides.  Clearly  front  equation  (1),  the  Lagrangian 
evolution  of  the  equatorial  zonal  flow  is  strongly  con¬ 
trolled  by  Xmech  variations  due  to  wave  forcing,  as 
models  of  the  quasi-biennial  oscillation,  SSAO,  and 
MLTSAO  have  demonstrated. 

However,  coherent  time  variations  in  wave  activity 
may  also  influence  U.  The  influence  that  diurnal  vari¬ 
ations  in  convectively  generated  gravity  waves  could 
have  on  tidal  structures  was  investigated  by  Eck- 
erniann  and  Marks  (1996).  They  showed  that  a  time 
variation  in  wave-field  intensity  yielded  a  time-varying 
Xmech  contribution,  whereupon  equation  (1)  can  be 
reexpressed  as 

DU  1  8PX 

n  —  XXtes,  (2) 

Dt  p  8t 

where  Xdiss  is  a  forcing  due  to  wave  dissipation  and 
turbulent  diffusion,  p  is  background  density  and  Px  is 
the  total  pseudomomentum  density  of  the  wave  field. 
Equation  (2)  clearly  shows  that,  even  when  small-scale 
dissipation  is  limited  (A,&1«0),  strong  Lagrangian 
accelerations  of  the  zonal  flow  can  result  from  local 
time  variations  in  the  total  wave  pseudomomentum 
density  impinging  upon  the  mesosphere. 

Here  we  associated  Px  with  the  total  pseudo- 
momentum  density  due  to  gravity  waves  and  diurnal 
tides  in  the  equatorial  MLT  region.  Our  analysis 
showed  strong  intraseasonal  variations  in  wave  vari¬ 
ances  and  tidal  amplitudes  (see  Figs  12  and  13).  Given 
realistic  assumptions  about  wave  field  anisotropy  and 
energy  distributions,  similar  periodicities  in  Px  may 
also  be  expected  (e.g.,  Fritts  and  VanZandt,  1993), 
whereupon  dPJdt  will  have  large  peaks  at  periods 
~  60  days,  ~ 40  days,  and  ~ 25  days.  From  equation 
(2),  we  see  that  for  the  simplified  situation  of  Xdiss  =  0 
and  D/Dtx8/8t ,  then  zonal  wind  accelerations  follow 
p~l8PJdt.  Thus  intraseasonsal  peaks  in  dPJdt  (e.g., 
Fig.  13)  can  appear  in  8U/8t  and  hence  in  U,  as 
observed  in  Fig.  4. 

Of  course  Xdlss=£0  and  D/Dt^=d/8t  at  MLT  heights, 
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so  that  while  transience  may  play  an  important  role, 
Xdiss  is  also  large  at  these  heights  and  will  play  a  pri¬ 
mary  role  in  the  mean-flow  evolution.  Dissipation  of 
intraseasonally  modulated  wave  and  tidal  activity  at 
MLT  heights  is  also  likely  to  produce  intraseasonal 
periodicities  in  Xdiss,  leading  via  equation  (2)  to  a  more 
complicated  mean-flow  evolution  (see  e.g.,  Section 
4.10  of  Baines  (1995))  which  should  also  result  in  an 
intraseasonal  variation  in  U,  as  argued  by  Eckermann 
and  Vincent  (1994).  Additional  influences,  such  as 
nonzonal  wave  breaking  and  meridional  transport 

^  jt  .  may  also  induce  similar  periodicities  in  v* 

to  those  observed  in  Fig.  4. 

While  we  have  presented  both  observational  and 
theoretical  evidence  to  support  gravity-wave  and  tidal 
driving  of  intraseasonal  oscillations  in  the  middle 
atmosphere,  more  detailed  numerical  modeling  will 
be  needed  to  assess  this  mechanism  fully.  Interestingly, 
Held  et  al.  (1993)  documented  a  60  day  stratospheric 
oscillation  that  arose  naturally  within  a  tropical 
numerical  model  which  contained  realistic  par- 
ameterizations  of  radiation,  moisture  and  convection. 
This  oscillation  was  presumed  to  be  forced  by  con- 
vectively  generated  gravity-wave  activity  in  the  model. 
While  these  model  results  and  our  observations  differ 
in  many  ways  (different  heights,  no  semiannual  oscil¬ 
lation  in  their  model,  and  so  on),  their  study  at  least 
establishes  a  precedent  for  the  forcing  of  intraseasonal 
wind  oscillations  in  the  tropical  middle  atmosphere 
by  convectively  generated  waves. 

Our  theoretical  discussion  has  assumed,  for 
simplicity,  that  the  peaks  in  Fig.  4  are  caused  by  mean 
wind  variations  over  intraseasonal  time  scales. 
Ground-based  data  such  as  these,  however,  cannot 
tell  us  whether  they  are  indeed  mean-flow  vacillations 
(wavenumber  0)  or  are  instead  due  to  planetary-wave 
disturbances  (wavenumbers  >  1).  In  fact,  either  is 
possible,  as  discussed  by  Eckermann  and  Vincent 
(1994).  Since  MJO-modulated  convection  occurs 
within  narrow  longitude  zones  rather  than  uniformly 
around  the  equator  (e.g.,  Dunkerton  and  Crum, 
1995),  this  will  yield  a  zonally  asymmetric  distribution 
of  gravity  waves  emanating  from  convective  zones. 
Simulations  by  Williams  and  Avery  (1996)  also  reveal 
longitudinal  variations  in  convectively  forced  tidal 
activity  in  the  equatorial  MLT  region.  This  zonally 
asymmetric  wave  activity  could  force  an  intraseasonal 
planetary  wave  in  the  MLT  region  (e.g.,  Holton, 
1984).  However,  similar  longitudinally  confined  con¬ 
vection  existed  in  the  model  of  Held  et  al.  (1993)  yet 
the  ensuing  60  day  stratospheric  oscillation  appeared 
in  the  mean  zonal  wind.  Since  equatorial  convection 


generates  waves  with  a  range  of  zonal  wavenumbers 
(e.g.,  Bergman  and  Salby,  1994),  then  the  mean-flow 
driving  produced  by  such  a  broad  spectrum  of  waves 
may  not  favor  the  generation  of  pure  planetary-wave 
modes  at  MLT  heights.  The  oblique  zonal  propa¬ 
gation  of  equatorial  modes  and  gravity  waves  away 
from  convective  zones  may  also  lead  to  a  more  sym¬ 
metric  zonal  distribution  of  wave  activity  at  heights 
some  distance  above  the  convection  (e.g.,  Alexander, 
1996)  and  thus  to  a  more  zonally  symmetric  flow 
forcing.  Since  there  are  also  strong  gravity  wave-tidal 
interactions  in  our  data  (see  also  Isler  and  Fritts,  1996; 
Connor  and  Avery  (1996)),  this  will  lead  to  mutual 
interactions  among  the  mean  flow,  the  intraseasonal 
flow,  tidal  oscillations  and  gravity  waves,  the  net 
results  of  which  are  difficult  to  predict  at  present. 
Satellite  observations  should  help  to  clarify  the  true 
zonal  structure  of  the  intraseasonal  MLT  winds. 


Origin  of  intraseasonal  periodicities  in  the  MLT  region 

One  difficulty  with  our  theory,  in  which  MJO- 
related  variations  are  forced  in  MLT  winds  through 
gravity-wave  and  tidal  coupling,  is  that  the  intra¬ 
seasonal  peaks  in  Fig.  4  differ  somewhat  from  those 
observed  in  the  lower  atmosphere.  The  MJO  peak  in 
the  troposphere  is  broad:  it  was  originally  named  the 
40-50  day  oscillation  by  Madden  and  Julian  (1971), 
although  later  studies  have  often  referred  to  it  as  the 
30-60  day  oscillation  to  stress  the  variable  and  broad¬ 
band  nature  of  the  peak  (see  the  review  of  Madden 
and  Julian  (1994)).  In  contrast,  the  peaks  observed  in 
Fig.  4  are  sharper,  and  are  displaced  somewhat  from 
the  mean  MJO  peak  at  ~45  days. 

The  band  B  peaks  at  35-40  days  seem  to  relate  most 
easily  with  the  30-60  day  MJO,  whereas  the  stronger 
quasi-60  day  peak  that  we  observe  in  band  A  is  at  the 
long-period  limit  of  MJO-related  variability. 
However,  the  activity  within  bands  A  and  B  seems  to 
occur  at  similar  times  (e.g.,  Fig.  5),  suggesting  that  the 
two  bands  are  closely  related.  Despite  considerable 
interannual  variability,  the  mean  seasonal  cycles  of 
activity  within  bands  A  and  B  showed  a  tendency  for 
an  annual  cycle,  with  greatest  activity  around  Jan- 
uary-February  and  least  activity  during  August  (Fig. 
7).  This  is  similar  to  both  the  interannual  variability 
and  mean  annual  variation  of  the  MJO  (see  e.g., 
Madden,  1986;  Hartmann  and  Gross,  1988;  Gutzler 
and  Madden,  1989),  and  is  consistent  with  the  hypo¬ 
thesized  MJO-modulated  tidal  and  gravity-wave  driv¬ 
ing  of  the  oscillations  within  bands  A  and  B. 

Some  clues  as  to  how  the  MJO  at  ~45  days  might 
drive  MLT  oscillations  of  ~  60  days  (band  A)  and 
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35M0  days  (band  B)  is  provided  by  the  periodicities 
evident  in  the  gravity-wave  variances  (Fig.  12).  At  86- 
90  km  there  is  a  zonal  peak  at  ~  60  days,  whereas  at 
92-98  km  the  zonal  peak  has  shifted  to  ~  45  days,  and 
peaks  in  meridional  velocity  variance  now  appear  at 
~  45  days,  ~  60  days,  and  ~  90  days.  We  postulated 
that  intraseasonal  variations  in  upper-level  meridional 
gravity-wave  activity  arose  on  interacting  with  the 
background  flow  below  90  km.  Below  90  km,  in 
addition  to  strong  quasi-60day  wind  variations  (of 
frequency  a>60).  there  is  of  course  a  strong  MLTSAO 
variation  of  ~  1 83  days  (frequency  eoSAO).  Thus,  grav¬ 
ity  waves  interacting  nonlinearly  with  these  two  flow 
periodicities  will  yield  a  time  series  of  gravity-wave 
variance  containing  not  only  the  frequencies  co60  and 
coSAO,  but  also  the  interactive  combinations  co60  +  coSAO 
These  two  combinations  give  rise  to  periods  of 
~  90  days  and  ~  45  days,  which  are  very  close  to  the 
peaks  observed  in  Fig.  12b. 

A  logical  extension  of  this  argument  may  explain 
how  the  band  A  and  B  peaks  are  produced  by  wave 
forcing  over  MJO-related  time  scales.  Since  the  MJO 
is  centered  at  a  period  of  ~  45  days  (frequency  coMJO), 
it  will  presumably  produce  a  similar  period  in  con- 
vectively  generated  gravity-wave  activity.  Gravity 
waves  are  believed  to  play  a  major  role  in  driving  the 
SSAO  (e.g.,  Sassi  et  al.,  1993;  Eckermann  et  al.,  1995). 
Thus,  when  this  MJO-modulated  gravity-wave 
activity  interacts  with  the  upper  stratospheric  flow  in 
forcing  the  SSAO,  the  emergent  gravity-wave  activity 
might  be  expected  to  contain  frequencies  coMJO  +  a>SAo- 
These  combinations  yield  periods  of  ~  60  days  and 
~  36  days,  which  are  very  similar  to  the  band  A  and 
band  B  peaks  detected  in  the  MLT  wind  oscillations 
(Fig.  4)  and  gravity-wave  variances  (Fig.  12a).  Similar 
effects  may  also  influence  diurnal  tidal  activity,  and 
may  explain  the  60  day  and  strong  40  day  peaks  in  this 
activity  observed  in  Fig.  13. 

The  observed  peaks  in  the  zonal  velocity  spectra 
within  the  20-25  day  range  (band  C)  seem  to  relate  to 
the  20-25  day  tropical  tropospheric  oscillation 
observed  in  the  western  Pacific  region.  This  oscillation 
is  known  to  induce  similar  periodicities  in  regional 
convective  activity  (Hartmann  et  al.,  1992;  Kiladis  et 
al.,  1994),  and  thus  (presumably)  in  the  intensity  of 
convectively  generated  gravity  waves  and  non¬ 
migrating  diurnal  tides  entering  the  middle  atmo¬ 
sphere.  Hartmann  et  al.  (1992)  noted  that  20-25  day 
TTIV  was  most  active  during  September-December, 
which  again  is  very  similar  to  the  seasonal  variations 
of  band  C  activity  observed  in  Fig.  7.  This  too  is  con¬ 
sistent  with  a  forcing  of  the  band  C  activity  by  con¬ 
vectively  generated  tides  and  gravity  waves.  Forcing 
by  the  diurnal  tide  seems  to  be  most  relevant  to  the 


band  C  activity,  since  there  is  a  strong  24  day  peak  in 
diurnal  tidal  activity  (see  Fig.  13). 


SUMMARY 

Our  analysis  of  over  five  years  of  horizontal  velocity 
measurements  of  the  MLT  region  over  Christmas 
Island  has  revealed  intraseasonal  variations  in  the 
zonal  winds,  as  initially  reported  by  Eckermann  and 
Vincent  (1994).  In  this  study,  we  have  also  found 
clear  intraseasonal  variations  in  the  meridional  winds, 
gravity-wave  activity,  and  diurnal  tidal  amplitudes, 
all  of  which  have  similar  periods  to  those  initially 
detected  in  the  zonal  winds. 

Three  main  period  bands  of  activity  in  the  zonal 
winds  were  identified,  at  around  60  days  (band  A), 
~ 40  days  (band  B),  and  ~  25  days  (band  C).  The 
60  day  signal  was  also  observed  strongly  in  gravity- 
wave  activity,  and  moderately  in  meridional  wind  and 
diurnal  tidal  activity.  A  35  day  signal  was  observed  in 
meridional  winds,  and  weak  40  day  and  25  day  peaks 
in  gravity-wave  activity  were  observed  at  86-90  km. 
However,  diurnal  tidal  amplitudes  exhibited  strong 
peaks  at  ~  40  days  and  ~  25  days  throughout  the  86- 
98  km  height  range.  Periodicities  in  gravity-wave 
activity  also  exhibited  height  variations  which  may  be 
due  to  nonlinear  coupling  with  the  background  flow. 

The  activity  in  all  three  intraseasonal  bands  was 
associated  with  intraseasonal  oscillations  of  the  trop¬ 
ical  troposphere.  The  peaks  within  bands  A  and  B 
were  within  the  range  of  the  30-60  day  Madden- 
Julian  Oscillation,  and  the  mean  annual  variation  of 
activity  in  bands  A  and  B  was  similar  to  that  observed 
for  the  MJO.  Band  C  activity  exhibited  a  period  range 
and  mean  seasonal  variation  that  agreed  well  with  20- 
25  day  oscillations  encountered  in  the  tropical  tropo¬ 
sphere. 

Our  observations  led  us  to  propose  that  30-60  day, 
20-25  day  and  diurnal  variations  in  tropical  tro¬ 
pospheric  convection  produce  intraseasonal  cycles  in 
the  intensity  of  gravity  waves  and  nonmigrating  diur¬ 
nal  tides  entering  the  middle  atmosphere.  Theoretical 
assessments  of  such  wave  production  indicated  that 
intraseasonal  variations  in  the  wave-induced  driving 
of  the  zonal  MLT  flow  will  result,  which  in  turn  will 
produce  intraseasonal  variations  in  the  mean  zonal 
winds,  as  observed.  Our  observation  of  similar  period¬ 
icities  in  diurnal  tidal  and  gravity-wave  activity  in  the 
MLT  region  support  this  theory. 

Further  observations  and  modeling  are  required  to 
understand  the  intraseasonal  behavior  of  the  equa¬ 
torial  MLT  region  better.  For  example,  it  is  still 
unclear  whether  these  intraseasonal  wind  disturbances 
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are  mean-flow  or  planetary- wave  phenomena.  It  is 
also  unclear  how  coherent  they  might  be  around  the 
equatorial  belt:  for  example,  we  showed  that  the 
activity  here  seems  to  differ  somewhat  from  intra- 
seasonal  winds  observed  near  the  stratopause  over 
India  (Ziemke  and  Stanford,  1991;  Kumar  and  Jain, 
1994).  Analysis  of  data  from  UARS  instruments  may 
provide  information  on  these  and  other  questions  that 


have  arisen  from  our  initial  analyses  of  equatorial 
MLT  wind  data  from  the  central  Pacific. 
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